Thermoelectric transport in Bi 2 Te3/Sb 2 Te3 superlattices 
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The thermoelectric transport properties of Bi2Te3/Sb2Te3 superlattices are analyzed on the basis 
of first-principles calculations and semi-classical Boltzmann theory. The anisotropy of the ther- 
moelectric transport under electron and hole-doping was studied in detail for different superlattice 
periods at changing temperature and charge carrier concentrations. A clear preference for thermo- 
electric transport under hole-doping, as well as for the in-plane transport direction was found for 
all superlattice periods. At hole-doping the electrical transport anisotropics remain bulk-like for all 
investigated systems, while under electron-doping quantum confinement leads to strong suppression 
of the cross-plane thermoelectric transport at several superlattice periods. In addition, insights on 
the Lorenz function, the electronic contribution to the thermal conductivity and the resulting figure 
of merit are given. 

PACS numbers: 31.15.A-,71.15.Mb,72.20.Pa,72.20.-i 
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I. INTRODUCTION 

Solid-state thermoelectric (TE) power generation de- 
vices possess the desirable nature of being highly reli- 
able, stable, compact and integrable and have potential 
applications in waste-heat recovery and outer space ex- 
plorations. However, while intensively studied in the last 
decades, poor energy conversion efficiencies below a few 
percent at room temperature prohibited the triumph of 
the TE devices as promising alternative energy sources. 
The conversion performance of a thermoelectric material 
is quantified by the figure of merit (FOM) 



ZT = 



(i) 



where a is the electrical conductivity, S the thcrmopower, 
n e i and K p h are the electronic and lattice contribution to 
the thermal conductivity, respectively. From Eq. 1 it is 
obvious, that a higher ZT is obtained by decreasing the 
denominator or by increasing the numerator, the latter 
being called power factor PF = aS 2 . While a, S, K e i 
and Kph can individually be tuned by several orders of 
magnitude, the interdependence between these properties 
impede high values for the FOM 1 ' 2 . Bi 2 Te 3 , Sb 2 Te 3 and 
their related alloys dominate the field of thermoelectrics 
with ZT around unity from the 1950's through today 3-5 . 

The idea of thermoelectric superlattices (SL) allows 
for concepts, which could enable both, the suppression 
of the cross-plane thermal conductivity 6-8 and the in- 
crease of the electronic power factor 9 . It suggests that 
cross-plane transport along the direction perpendicular 
to the artificial interfaces of the SL reduces phonon heat 



conduction while maintaining or even enhancing the elec- 
tron transport 10 . In 2001 a break-trough experiment 
by Venkatasubramanian et al. reported a record appar- 
ent ZT = 2.4 for p-typc Bi 2 Te3/Sb 2 Te3 and ZT = 1.4 
for n-type Bi 2 Tc 3 /Bi 2 Te 2 .83Seo.i7 superlattices 8 ' 11,12 , al- 
though this values have not yet been reproduced to the 
best of our knowledge. 

With the availability of materials with ZT > 3 ther- 
moelectric materials could compete with conventional 
energy conversion methods and new applications could 
emerge 13 . Beside thermal conductivities below the al- 
loy limit, the investigations of Venkatasubramanian et 
al. 12 found a strong decrease of the mobility anisotropy 
and the related electronic thermoelectric properties for 
the SLs at certain periods. This is counter-intuitively, as 
superlattices are anisotropic by definition and even the 
tclluridc bulk materials show intrinsic anisotropic struc- 
tural and electronic properties 14-18 . While considerable 
effort was done in experimental research 19-26 , theoretical 
investigations on Bi 2 Te3/Sb 2 Te3 SLs are rare. Various 
available theoretical works concentrate on the electronic 
structure and transport properties of the bulk materi- 
als 17.27-39 with 

some of them discussing the influence 
of strain, which could occur at the SL interfaces 18,30 ' 31 . 
To our knowledge, only a sole theoretical work discussed 
the possible transport across such tclluridc SL structures. 
Based on density functional theory, Li et al. 32 focussed 
on the calculation of the electronic structure for two dis- 
tinct Bi 2 Te3/Sb 2 Te3 SLs, stating changes of the mobility 
anisotropy estimated from effective masses. To extend 
this work and to clarify the open questions on the re- 
duced mobility anisotropy and the enhanced thermoelec- 
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trie efficiency, we are going to discuss in this paper the 
anisotropic thermoelectric electronic transport of seven 
different Bi 2 Te 3 /Sb2Te3 SLs, including the bulk materi- 
als, on the basis of density functional theory and semi- 
classical transport calculations. 

For this purpose the paper will be organized as fol- 
lows. In Sec. II we introduce our first principle electronic 
structure calculations based on density functional theory 
and the semi-classical transport calculations based on the 
solution of the linearized Boltzmann equation. A brief 
discussion of the obtained band structures, including the 
dependence of the band gap on different SL periods, is 
done in section Sec. III. With this, we present in Sec. IV 
the electronic thermoelectric transport properties, that 
is electrical conductivity, thermopower and the related 
power factor, of the Bi 2 Te3/Sb2Te3 SLs at different SL 
periods with a focus on their directional anisotropics. 
The discussions cover a broad temperature and doping 
range and will conclude the bulk materials. Even though 
only p-type conduction was found in SLs stacked out 
of pure Bi2Te3 and Sb2Te3 n-type conduction will be 
studied, too, as being possible at appropriate extrinsic 
doping. For the latter case strong quantum confinement 
effects were found, which will be discussed in detail in 
Sec. IV B. To give a clue on possible values for the figure 
of merit, in Sec. IV C results for the electronic contribu- 
tion to the thermal conductivity, the Lorenz function as 
well as existing experimental results for the lattice part 
of the thermal conductivity will be presented. Most of 
the discussions are done in a comparative manner con- 
sidering the experimental findings of Rcf. 12 . 



II. METHODOLOGY 

For both bismuth and antimony tclluridc, as well as for 
the composed heterostructures, we used the experimen- 
tal lattice parameters and relaxed atomic positions 33 as 
provided for the hexagonal Bi2Tc3 crystal structure with 
15 atomic layers. The layered structure itself is com- 
posed out of three formula units, Tci-Bi-Tc2-Bi-Tei, of- 
ten called quintuples. The hexagonal lattice parameters 
are chosen to be a|f£ e = 4.384A and c^ffte = 30.487A 
for Bi 2 Te 3 , Sb 2 Te 3 and the BiaTeg/SbaTeg SLs, respec- 
tively. In fact, the main difference between the lattices 
of Bi 2 Te3 and Sb2Te3 is a decrease of the in-plane lattice 
constant with an accompanied decrease in cell volume. 
So, a change between the two lattice constants can be 
related to either compressive or tensile in-plane strain. 
Preceding intense studies revealed that a larger in-planc 
lattice constant, e.g. Og?^ e > agj^. e , is favourable for an 
enhanced cross-plane TE transport 18,29,34 . For this pur- 
pose, the experimental lattice parameters of Bi 2 Te3 were 
chosen for the studied heterostructures. Structural re- 
laxations revealed only minor influences on the bulk 
electronic structure 29-32 and are beyond the scope of 
this work. To introduce SLs with different layer peri- 




FIG. 1. (color online) Shown are three possible hexagonal 
unit cells of the (Bi2Te3) I /(Sb2Te3)i_ 2! superlattices. (a) x=l 
which is bulk Bi2Te3, (b) x= 2 /3 and (c) x= x /6. A concentra- 
tion x=0 would coincide with bulk Sb2Te3 in the lattice of 
Bi 2 Te 3 . 



ods compared to the experiments of Vcnkatasubrama- 
nian et al. 12 we subsequently substitute the Bi site by Sb, 
starting with six Bi sites in hexagonal bulk Bi2Te3 (see 
Fig. 1(a)). For instance, substituting two atomic layers 
of Bi with Sb leads to a {Bi 2 Te 3 ) x / (Sb 2 Te 3 ) 1 _ x SL with 
x = 2 /3, that is two quintuple Bi2Te3 and one quintuple 
Sb2Te3 (see Fig. 1(b)). The latter case coincides with 
a (20A/10A)-(Bi 2 Te3/Sb2Te 3 ) superlattice in the exper- 
imental notation of Ref. 12 . 

Our thermoelectric transport calculations are per- 
formed in two steps. In a first step, the detailed elec- 
tronic structure of the Bi2Te3/Sb2Te3 SLs were ob- 
tained by first principles density functional theory cal- 
culations (DFT), as implemented in the fully relativis- 
tic screened Korringa-Kolm-Rostoker Greens-function 
method (KKR) 35 . Within this approach the DiRAC- 
cquation is solved sclf-consistently and with that spin- 
orbit-coupling (SOC) is included. Exchange and cor- 
relation effects were accounted for by the local density 
approximation (LDA) parametrized by Vosco, Wilk, and 
Nusair 36 . Detailed studies on the electronic structure and 
transport anisotropy of the bulk tellurides Bi2Te3 and 
Sb2Te3 were published before 18,34 and show very good 
agreement to experimental results and other theoretical 
findings. 

With the well converged results from the first step we 
obtain the thermoelectric transport properties by solving 
the linearized Boltzmann equation in relaxation time ap- 
proximation (RTA) within an in-house developed Boltz- 
mann transport code 3 ,38 . Boltzmann transport calcula- 
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tions for thermoclcctrics have been carried out for quite a 
long time and show reliable results for wide- and narrow 
gap semiconductors 38-42 . Calculations on the electronic 
structure and TE transport for bulk Bi 2 Te 3 17 ' 31 ' 43 ' 44 and 
Sb2Te3 28,31 ' 45 were presented before. Here the relaxation 
time r is assumed to be isotropic and constant with re- 
spect to wave vector k and energy on the scale of ksT. 
This assumption is widely accepted for degenerate doped 
semiconductors. Within the RTA, from comparison of 
the calculated electrical and electronic thermal conduc- 
tivities (eq. 3 and 4) with experiment it is possible to 
conclude on the relaxation time. In the following r is set 
to 10 fs for all bulk and heterostructure systems, regard- 
less of any directional anisotropy or charge carrier depen- 
dence. We note, that within RTA for the thermopower 
S (eq. 3) the dependence of the transport distribution 
function (TDF) , as introduced in the next paragraph, on 
the energy is essential. That is, not only the slope of 
the TDF, moreover the overall functional behaviour of 
the TDF on the considered energy scale has to change to 
observe an impact on the thermopower. 

Within the RTA the TDF (fx, 0) 46 and with this 

the generalized conductance moments £}-"',,(fj,,T) are de- 
fined as 



(2 



^(^,(x,||)) 2 (^-M)"(- 



dE 



E=EY 



(2) 

E% denotes the band structure of band v, v% the group ve- 
locity and f{n,T) the FERMi-DiRAC-distribution. v v k 
v k (_L) denote the group velocities in the directions in the 
hexagonal basal plane and perpendicular to it, respec- 
tively. Within here the group velocities were obtained 
as derivatives along the lines of the Blochl mesh in the 
whole Brillouin zone (BZ) 34 . The directions of these lines 
are parallel to the reciprocal space vectors and so the 
anisotropy of the real lattice is reflected in these vec- 
tors. A detailed discussion on implications and difficul- 
ties on the numerical determination of the group veloc- 
ities in highly anisotropic materials was currently pub- 
lished elsewhere 47 . As can be seen straight forwardly, 
the temperature- and doping-dependent electrical con- 
ductivity a and thermopower S in the in- and cross-plane 
directions arc defined as 



e 2 C 



(0) 

-J--II 



s 



1 gSiO^T) 



(3) 



and the electronic part to the total thermal conductivity 
accounts to 



«eJ_L,|| 



^(4 2 ,ii(^ T ) 



4°ji(M.T) 



(4) 



The second term in eq. 4 introduces corrections due to the 
Peltier heat flow that can occur when bipolar conduction 
takes place 48 . 



The chemical potential fj, at temperature T and extrin- 
sic carrier concentration N is determined by an integra- 
tion over the density of states (DOS) n(E) 



VBM 



N = 



J &En{E)[f{n,T)-l]+ J dEn(E).f(p,T$) 



CBM 



where CBM is the conduction band minimum and VBM 
is the valence band maximum. The k-spacc integration 
of eq. 2 for a system with an intrinsic anisotropic texture, 
e.g. in rhombohedral and hexagonal structures, is quite 
challenging. In preceding publications 34 ' 47 we stated on 
the relevance of adaptive integration methods needed to 
reach convergence of the energy dependent TDF. Espe- 
cially in regions close to the band edges, which are evi- 
dent for transport, the anisotropy of the TDF requires a 
high density of the k-mesh. Here, convergence tests for 
the transport properties showed that at least 150 000 k- 
points in the entire BZ had to be included for sufficient 
high doping rates (N > 1 x 10 19 cm -3 ), while for energies 
near the band edges even more than 56 million k-points 
were required to reach the analytical effective mass val- 
ues and the corresponding conductivity anisotropies at 
the band edges. 



III. ELECTRONIC STRUCTURE 

In Figs. l(a)-(g) the electronic bandstructures 
on the hexagonal high symmetry lines for all 
(Bi2Te3)a;/(Sb2Te3)i_ x SLs are shown, starting with (a) 
x = which is tensile strained bulk Sb2Te3 and ending 
with (g) x = 1, which is bulk Bi 2 Te3 . For the case of 
(b) x = 1 /6, (d) x = 1 /2 and (f) x — 5 /6 a further band 
splitting can be noticed, which stems from the missing 
space inversion symmetry in these systems and with that 
the former band degeneracy is lifted. This situation al- 
ways occurs if the Bi(Sb) sites in each quintuple are not 
uniformly occupied. 

For increasing number of Bi layers in the SLs no dras- 
tic change in the band structure topology can be stated. 
Only slight variations were found for the in-plane band 
directions. Of stronger impact could be the change of 
band dispersion which occurs for the lowest lying con- 
duction band in the cross-plane direction TA. Here a con- 
tinuous change of the bands slope is found for increasing 
amount of Bi layers in the SL. An almost vanishing dis- 
persion and very flat bands in cross-plane direction are 
found for the SL with x = 0.5, which is three Bi-like lay- 
ers and three Sb-like layers. Further amount of Bi layers 
in the system leads to an increase in the bands slope, 
while showing different sign compared to bulk Sb2Te3. 

In Fig. 2(h) the calculated band gap in dependence 
on the superlattice period is shown. Applying an ex- 
tended tetrahedron method 50 ' 51 and very dense k-mesh's 
in the BZ, the band gap values were determined within 
an uncertainty below 1%. While for Sb2Te3 at the ex- 
perimental lattice parameters, we previously found a di- 
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FIG. 2. (color online) Electronic bandstructures for (Bi2Te3) ;E /(Sb2Te3)i_ ;E superlattices in the hexagonal unit cell with different 
superlattice periods, (a) x = 0, (b) x = (c) x = 1 /z, (d) x = 1 /2, (e) x — 2 /s, (f) x = 5 /6 and (g) x = 1. In (h) black 
circles show the calculated fundamental energy gap in dependence on the superlattice period. Red squares show experimental 
findings 49 for (Bi2Te3) I /(Sb2Te3)i_ a; alloys, which were linearised to allow for comparison with our LDA results. 



rect band gap located at the center of the BZ 3 , an 
indirect gap of E g = 140 mcV can be stated for the in- 
plane tensile strained Sb2Te3. For Bi2Tc3 as well as for 
all Bi2Tc3/Sb2Te3 SLs indirect band gaps are obtained, 
too. A known difficulty within standard DFT is the gen- 
eral underestimation of the semi-conductors band gaps 
at zero temperature 52 , as well as the missing tempera- 
ture dependence of E g (T) 53 . For small band gap ther- 
moelectrics, such as Bi2Te3 and Sb2Te3, this could im- 
pinge the TE transport. The thermopower might be re- 
duced at high temperature and low doping due to bipo- 
lar conduction 17,39 . This effect would be overestimated 
if the band gap is underestimated. With E g (T) being 
considered for TE bulk materials 17 ' 38 lack of knowledge 
on the absolute size of the gap, as well as its temperature 
dependence permits such gap corrections for the strained 



bulk materials as well as for the SLs. However, the calcu- 
lated bulk band gap of E g — 105 meV for unstrained bulk 
Bi2Te3 is in better agreement with the experimental value 
of E g = 130 me V 49 , than for unstrained Sb2Te3 where a 
calculated value of E g = 90 mcV faces experimental val- 
ues between E g = 150 — 230 me V 49,54 . As is well known, 
Bi2Te3 and Sb2Te3 exhibit band inversions at certain ar- 
eas in the BZ 55 . Within LDA the strength of band in- 
version is most likely underestimated 56 . At a given band 
inversion the strength of the spin orbit interaction then 
controls the size of the band gap. Fortuitously, overes- 
timated SOC effects and underestimated band inversion 
tend to cancel each other leading to good results for the 
band gap size and wave function character. 

As a lack of data, we can only compare the 
Bi 2 Te 3 /Sb2Te3 SL with x = 1/3 (see Fig. 2(c)) with pre- 
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vious results of Li et al. 32 . While they applied a full- 
potential Linearized Augmented Plane Wave (FLAPW) 
method and treated spin-orbit coupling as a second or- 
der perturbation, the results on the bands topology are in 
very good agreement. However, our band gap is substan- 
tially larger with E g = 95 meV neglecting any structural 
relaxations, compared to their value of E g = 27meV. 



IV. THERMOELECTRIC TRANSPORT 

A. Effects of superlattice period composition 

With knowledge on the electronic structure, we are 
now able to calculate the desired thermoelectric trans- 
port properties of the Bi 2 Te 3 /Sb2Te3 SLs. As a starting 
point the in-plane transport properties at room temper- 
ature for the electron doped (cf. Figs. 3(a)-(c)) and hole 
doped (cf. Figs. 3(d)-(f)) hctcrostructures are chosen, 
while afterwards the anisotropy referring to the transport 
in cross-plane direction is discussed in detail. Prelimi- 
nary studies revealed the optimal charge carrier concen- 
tration for the SLs to be in the range of 3 — 6 x 10 19 cm -3 
18 . For the sake of clarity results are presented for 
three relevant charge carrier concentrations of 3, 6 and 
9 x 10 19 cm~ 3 (cf. solid, dashed and dotted lines in Fig. 3, 
respectively). 

Under electron doping (cf. Fig. 3(a)) a decrease of 
the in-plane electrical conductivity for the superlattices 
compared to the bulk materials is found. This decrease is 
more pronounced at higher charge carrier concentrations, 
while only slightly being dependent on the SL period. 
Despite taking into account an isotropic and constant re- 
laxation time r = 10 fs, we find very good agreement with 
experiment 57 ^ 58 for bulk Bi 2 Te 3 with a = 1030 (fkm) -1 
at N = 3 x 10 19 cuT 3 . 

The absolute value of the n-type in-plane thermopower 
is shown in Fig. 3(b). At a carrier concentration of 
N = 3 x 10 19 cm~ 3 a higher amount of Bi2Te3 in the 
superlattices leads to a monotonically increase in the 
thermopower from Sii = 103/iV/K (bulk Sb 2 Te 3 ) to 
5|| = 141 //V/K (bulk Bi 2 Te3), while showing a dip at 
a composition of x = 2 /6 with S\\ below 100/iV/K. The 
latter anomaly is linked to confinement effects and is dis- 
cussed in detail in Sec. IV B. This overall behaviour of 
S,, is retained for higher charge carrier concentrations 
at reasonable smaller absolute values. Assembling the 
previous results, the power factor PF,, under relevant 
electron doping is shown in Fig. 3(c). Clearly, the re- 
duction of in-plane electrical conductivity cr. , as well as 
the dip of S,, at a SL period of x — 2 /& lead to a mini- 
mal power factor of about 4 /iW/cmK 2 at the named SL 
period. We find PF.. for the SLs always to be smaller 
than expected from an interpolation of the bulk abso- 
lute values. Furthermore due to compensation effects 
of and the dependence on the amount of dop- 
ing is less drastically for PF,, than for it's constituents. 



The best power factor was found for bulk Bi 2 Tc3 to be 
PF,, = 21 /iW/cmK 2 , while experimentally thin films 
and single crystals show PF y = 8 — 27//W/cmK 2 and 
PF = 45/xW/cmK 2 , respectively 25,59 . We note here, 
that in experiment n-type conduction was only apparent 
for Bi 2 Te3/Bi 2 Te 2 .83Seo.i7-SL. Nevertheless, to get more 
insight the physical mechanisms in thermoelectric SL 
transport, n-type transport in Bi 2 Te3/Sb 2 Te3 SLs should 
be of enhanced interest, too. 

Highest power factors and FOM were experimentally 
found for p-type Bi 2 Te3/Sb 2 Te3 SLs. The preference for 
hole conduction is dedicated to the large inherent defects 
introduced by the Sb 2 Te3 layers. In Figs. 3(d)-(f) the 
in-plane thermoelectric transport properties under hole 
doping are displayed in the same manner as done be- 
fore. Compared to the electron doped case (cf. Fig. 3(a)) 
the hole electrical conductivity cr (| is higher at the same 
charge carrier concentration. Furthermore almost no de- 
crease of cr (| could be found for the SLs, while this is more 
visible at lower charge carrier concentrations. For the 
in-plane thermopower the values at different superlattice 
compositions are again only slightly suppressed compared 
to the bulk systems. For a (Bi 2 Tc3) 2; /(Sb 2 Te3)i_ 2: SL at 
x = 3/ 6 we state S\\ = 149yuV/K, while S\\ = 154/iV/K 
and 5|| = 150/iV/K were found for bulk Sb 2 Tc3 and 
Bi 2 Tc3 at the lowest charge carrier concentration, re- 
spectively This negative bending of the thermopower at 
different superlattice periods is reflected and enhanced 
for PF y . From Fig. 3(f) it can be seen, that the in-plane 
power factor PF,, for the various superlattice is decreased 
compared to the bulk materials. However, the largest 
suppression (x = 2 /& and x = 3 /6) is found to be about 
20% compared to the bulk values, but still offers thermo- 
electric feasible values about PF,, = SO/xW/cmK 2 . 

To give a reference, in Table I the calculated in-plane 
thermoelectric properties are compared to experimental 
results. In the original work of Vcnkatasubramanian et 
al. 12 very large values of cr,, and result in a huge power 
factor PF|| about 72 yuW/cmK 2 at room temperature . 
These reported values are way larger than found for bulk 
or thin film Bi 2 Te 3 , Sb 2 Tc3 or their related alloys 25 . 
However, in a more recent study Winkler et al. 26 reported 
values for the in-plane electrical conductivity and ther- 
mopower of a comparable (Bi . 2 Sb .8) 2 Te3/Sb 2 Te3 sput- 
tered SL (cf. Table I), which are in very good agreement 
to our theoretical calculations and combine to an in-plane 
power factor PF,, above 30/xW/cmK 2 . This is similar to 
values for bulk single crystals with comparable composi- 
tions. In contrast to the original experiments 12,61 , which 
used low-temperature metal-organic chemical vapor de- 
position (MOCVD), Winkler et al. applied the concept 
of "nano- alloying" 21 . Here the elemental layers Bi, Sb, 
and Te are deposited by sputtering and subsequently an- 
nealed to induce interdiffusion and a solid-state reaction 
to form the SLs. The pronounced periodicity and c- 
orientation of the SLs have been demonstrated by sec- 
ondary ion mass spectrometry (SIMS) and X-ray diffrac- 



6 



g 3000 
O 

2000 



(a) 


1 , | 




3 . 10 8 

6 10 » -e- 
9-10 18 •• •»•• 




.A-'*' J 


'-A- 


A-"" 


A ^ 

Q * 


...A-"' 
Q-~~~ 


.J 

-0-''' 


r — -Q- 


a • 








: 


■ 


i 


■ 





20 



5 15 



0.2 0.4 0.6 O.i 

( Bi 2 Te 3)x / (Sb 2 Te 3 ) 1 . x 












1 


-(c) 








/ \ 


















- - - -A" * 




■ 


■ 


■ 


■ 





O 3000 
S 

o 

rt 2000 



(d) 



| 

-3. 100 



C/3 



0.2 0.4 0.6 

(Bi 2 Te 3 ) x /(Sb 2 Te 3 ) 1 



1 


■ 






1 






Q___ 


e 




D 1 

-9 ' 


■ 

■(e) 


A- - - - 

■ 


A 

■ 


A- 

■ 


-A- ' 




FIG. 3. (color online) Absolute values of in-plane thermoelectric transport properties for (Bi2Te3) a; /(Sb2Te3)i_ :E superlattices 
in dependence on the superlattice period. Shown are (a),(d) electrical conductivity (b),(e) thermopower Sy and (c),(f) power 
factor PF|| . The temperature is fixed to 300 K and results for three different charge carrier concentrations (in units of cm" 3 ) 
are stated, (a), (b), (c) refer to electron doping, while (d), (e), (f) refer to hole doping. 



TABLE I. Theoretical and experimental in-plane thermoelec- 
tric properties of p- type (Bi2Te3) :E /(Sb2Te3)i_ a ; SLs at room 
temperature. The materials composition amounts in all con- 
sidered systems to about x = See text for additional 
details. 



N 


a 


5|| 


PP„ 


Ref. 


(10 19 cm- 3 ) 


(Ocm)- 1 


GuV/K) 


(pW/cmK 2 ) 




3.0 


1300 


151 


30 


this work 


3.1 


1818 


~ 200 


~ 72 


[12, 61-63] 


3.2 


761 - 1160 


172 - 189 


27-34 


[26] a 


5.8 


3050 


115 


40 


[64] b 



a sputtered (Bi .2Sb .8)2Te 3 /Sb2Te3 SL 
b (Bi^Sbi-a^Tes mixed crystal 

tion (XRD), respectively. 

While up to now we considered only in-plane trans- 
port, in the following the cross-plane transport of the 
superlattices will be discussed. The transport direction 
is therefore along the SL direction, perpendicular to the 
hexagonal basal plane of the bulk materials. In detail the 
directional anisotropy of the transport properties at room 
temperature are depicted in Figs. 4(a)-(c) and (d)-(f), for 
electron and hole doping, respectively. To get the abso- 
lute values for cross-plane transport, the in-plane val- 
ues previously shown in Fig. 3 should be divided by the 



anisotropics presented hereinafter. Anisotropics larger 
than unity represent suppressed thermoelectric transport 
in cross-plane direction and are therefore less desirable. 
As has been previously proven by experiment 14-16 ' 44,57 
and theory 17,18 ' 27,28 , already the bulk thcrmoelectrics 
Bi2Te3 and Sb2Te3 show large anisotropies for the elec- 
trical conductivity, thermopower and the related power 
factor. 

For a sense of purpose, the thermoelectric transport 
anisotropies under influence of hole doping will be con- 
sidered first. In Fig. 4(d) the anisotropy ratio of the elec- 
trical conductivity for various SL periods is illustrated 
at a temperature of 300 K. The anisotropy a \\/a 1 _ de- 
velops smoothly and monotonously between the bulk 
limits of a \\ja 1 _ = 2.7 and about < 7 i\/a ± = 5 — 6 for 
bulk Sb2Te3 and Bi 2 Te3, respectively. With increasing 
amount of Bi 2 Te3 in the superlattices the dependence of 
a ii/(7± on the charge carrier concentration is more pro- 
nounced. This is in accordance to previous findings for 
the bulk materials 18 . For the thermopower anisotropy 
s ii/s_l this picture holds, too. While for Sb2Te3 only a 
slight anisotropy of about s \\/s± =0.9 is found, the asym- 
metry increases for increasing amount of Bi in the SLs, 
saturating to about s i\/s± =0.75 for bulk Bi 2 Te3. The 
fact of the cross-plane thermopower being enhanced com- 
pared to the in-plane part is well known for the two bulk 
tellurides and compensates somewhat the high electri- 
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FIG. 4. (color online) Directional anisotropies of thermoelectric transport properties for (Bi2Te3) x /(Sb2Te3)i_ :E superlattices 
in dependence on the superlattice period. Shown are (a),(d) electrical conductivity ratio tT \\/rr± (b),(e) thermopower ratio 
S ll/Sj_ and (c),(f) power factor ratio pf ii/pf_,_. The temperature is fixed to 300K and results for three different charge carrier 
concentrations (in units of cm -3 ) are compared, (a), (b), (c) refer to electron doping, while (d), (e) and (f) refer to hole doping. 



cal conductivity anisotropy a \\/<j± to result in a less sup- 
pressed cross-plane power factor 18,44 . The anisotropy for 
the latter is shown in Fig. 4(f). Obviously, pf ii/pf_l is 
well above unity for all systems indicating a less pre- 
ferred cross-plane electronic transport. Compared to the 
bulk values of pF||/pf_l = 2.2 and pF||/pf_l = 2 - 3.5 for 
bulk Sb2Te3 and Bi2Te3 , respectively, the power factor 
anisotropy is only slightly larger for the SLs with differ- 
ent periods. As an example, for x = ^/e, which refers 
to a 5A/25A Bi 2 Te 3 /Sb2Te3 SL, the cross-plane power 
flow pf ii/pf_l is only suppressed by 13% with respect to 
bulk Sb2Tc3, while being enhanced by 26% compared to 
bulk Bi2Te3 at the optimal charge carrier concentration 
of N = 3 x 10 19 cm~ 3 . For the thermal conductivity in 
the SLs a suppression compared to bulk and the related 
alloys by about a factor of five is expected 8,12 . This would 
clearly lead to a benefit for the resulting FOM in compari- 
son to bulk, which is discussed more in detail in sec. IV C. 
However, we want to mention, that in the experiments of 
Vcnkatasubramanian et al. a further decrease in the elec- 
trical conductivity anisotropy was found for thin SLs at 
various SL periods 12 . It was stated that a \i/o-± under hole 
doping is about 0.8 - 1.4 for the Bi 2 Te3/Sb 2 Te3 SLs at 
different SL periods and therefore electrical cross-plane 
transport is strongly improved compared to bulk. Our 
calculations do not show such a trend. 



"\\/(t± under electron doping is shown. Bulk Bi2Te3 and 
bulk Sb2Te3 show anisotropies around a \\/a± = 8 and 
a \\/ax = 5, respectively, in good agreement to earlier 
studies 18 . One easily recognizes the anisotropy ratios to 
be larger than expected from the two bulk limits, while 
obtaining substantially large values of about a u/a ± = 20 
at TV = 3 x 10 19 cm~ 3 for a SL period of x = 
that is a 10A/20A Bi 2 Te 3 /Sb2Te3 SL, or one quintu- 
ple of Bi 2 Tc 3 and two quintuples of Sb 2 Te 3 . For the 
considered case the anisotropy strongly depends on the 
amount of doping, while decreasing rapidly at increased 
charge carrier concentration, but still reaching a \\/<j ± > 10 
at N — 9 x 10 19 cm~ 3 . At the same time the ther- 
mopower anisotropy shows a clear cross-plane preference 
at s \i/s± « 0.5. Nevertheless, the resulting power fac- 
tor anisotropy shows disappointing high values of about 
pf ii/pFj_ 3> 5 for the distinct SLs, while even showing 
PF \\/pF ± about 20 for the SL at a composition of x = 
2 /6. The suppressed cross-plane thermoelectric transport 
can clearly be linked to the large electrical conductivity 
anisotropics found for the n-type Bi2Tc3/Sb2Te3 SLs. In 
the following we want to emphasize, that the latter are 
related to quantum well effects in the conduction band, 
which are evoked by a conduction band offset between 
Bi 2 Te 3 and Sb 2 Te 3 in the SLs. 



In Fig. 4(a) the electrical conductivity anisotropy 
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B. Quantum well effects 

In the early 1990's concepts were presented to en- 
hance in-plane thermoelectric properties due to the use 
of quantum-confinement effects in SLs 65,66 . While huge 
enhancements on the in-plane figure of merit were pre- 
dicted, the authors suppressed electron tunnelling and 
thermal currents between the layers by introducing in- 
finite potential barriers and zero barrier widths. Later 
on it was shown, that for realistic barrier heights and 
widths the enhancement is rather moderate, predicting 
ZT values that at theirs best are a few percent larger 
than corresponding bulk materials 67,68 . 

It is known for the two tellurides, that due to SOC 
induced band inversions near the T point the CBM is 
derived from states localized at the Te atoms whereas 
the VBM is formed by Bi or Sb orbitals 55,69 . Neverthe- 
less, contributions from other areas in the BZ lead to 
the fact, that the valence band conduction contribution 
is mainly determined by Te states at appropriate charge 
concentrations. Due to this, a substitution of Bi and 
Sb in Bi2Te3 or Sb2Te3 affects the character of the va- 
lence band states only marginally and almost bulk-like 
electronic structure and transport properties can be ex- 
pected in the SLs. 

Venkatasubramanian et al. argued that in the SLs 
due to weak-confinement and near-zero band-offset, there 
is minimal anisotropy between in-plane and cross-plane 
electrical conductivities. While we can confirm that the 
valence band offset is almost vanishing in all SLs, the 
argumentation of Ref. 12 would conclude, that in bulk 
Bi2Te3 and Sb2Te3 < T \\/<r± = 1, as the band offsets in bulk 
materials are zero by definition 70 . The latter conclusion 
is obviously not the case. 

Bi2Te3 and Sb2Te3 show a theoretical band gap dif- 
ference of about 35 meV, thus a band offset in the SL is 
expected. Our calculations reveal that this difference is 
mostly located in the conduction bands. The offset in 
the conduction band edges of Bi2Te3 and Sb2Te3 sets 
up potential barriers in the superlattice, which leads 
to confinement of the electrons in the well regions. 
Therefore, the site resolved probability amplitude for 
the two bulk tellurides, as well as for the superlattice 
which showed the highest conductivity anisotropy, i.e. 
10A/20A Bi 2 Te 3 /Sb2Te3 SL (x = 2/e), is shown in Fig- 
ure 5(a)-(c). 

For Bi 2 Te3 the conduction band edge is dominated by 
contributions of the Tei followed by the Bi sites. A local- 
ization of states nearby the van der Waals gap is al- 
ready visible in the bulk system. The fact, that the band 
inversion does not heavily affect the orbital character at 
the CBM is caused by the indirect band gap character of 
Bi2Te3. Areas of the BZ where no band inversion occurs 
dominate the CBM. 

For Sb2Te3 this is quite different. Even though 
Sb2Te3 changes from an direct to an indirect semicon- 
ductor under applied in-plane tensile strain 34 , the CBM 
remains nearby the T point in the BZ. Therefore most 




FIG. 5. (color online) Site resolved probability ampli- 
tude for (a) bulk Bi 2 Te 3 , (b) bulk Sb 2 Te 3 and (c) a 
(Bi2Te3) a; /(Sb2Te3)i_ a ; superlattice with x = 2 /e. In the 
same manner the site resolved cross-plane electrical conduc- 
tivity a ± (in units of (Slcm) -1 ) is shown for (d) bulk Bi2Te3 , 
(e) bulk Sb 2 Te 3 and (f) a (Bi 2 Te 3 ) :c /(Sb2Te3)i-x superlat- 
tice with x = 2 /e. The temperature is fixed to 300 K and the 
charge carrier concentration is set to N — 3 x 10 19 cm -3 . 



of the contributions to conductivity arise from an area 
around the T point, where the band inversion favours the 
Sb character. As can be seen from Figure 5(b) the contri- 
butions to the DOS are equally distributed over all posi- 
tions in the unit cell, while slightly being enlarged on the 
Sb sites. For the (Bi 2 Te 3 ) a; /(Sb2Te3)i- K SL at x = 2 /e 
(cf. Figure 5(c)) we found quantum well states, which 
localize about half of the density in the Bi 2 Tc3 quintu- 
ple, while the density in the two Sb2Te3 quintuples is 
strongly depleted. We obtain similar results, if two quin- 
tuples (x — 4 /6) arc occupied by Bi2Te3. 

As shown in Figure 5(f) this quantum confinement is 
reflected in the contribution to the cross-plane electri- 
cal conductivity. Here, the local cross-plane conductivity 
c_i_(r) is calculated as introduced in Ref. 71, by weight- 
ing the contributions to C±u (m> ^) with the normalized 
probability amplitude ^(r)! 2 of the electronic states at 
chemical potential [i. Summing up o±(r) over all sites 
gives the total electrical conductivity <r± (cf. Fig. 3(a) 
and Fig. 4(a)). Weighting the DOS n(//) with the nor- 
malized probability amplitude refers to the local 
DOS n(p t r) (LDOS). 

Nevertheless sites in the Bi2Te3 quintuple with more 
accumulated density carry a larger contribution to the 
conductivity a ± , the total contribution compared to the 
bulk tellurides is strongly suppressed (comp. Figure 
5(d),(e)). This can be affirmed by a picture that elec- 
trons travelling in the cross-plane direction are exposed 
to a tunneling-like behaviour for about the distance of 
the Sb 2 Tc3 quintuples. This clearly leads to a diminished 
cross-plane group velocity of the electronic states. Com- 
paring Figure 5(d) and (e) we see furthermore that even 
the localization inside the quintuple in bulk Bi2Te3 can 
lead to reduced cross-plane electronic transport, reflected 
in larger total anisotropies about cr \\/cr ± = 8 for bulk 
Bi 2 Te3 compared to < T \\/<r± = 5 for bulk Sb2Te3. 

To extend to the results obtained at room temper- 
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ature in Figure 6 the temperature dependence of the 
in-plane and cross-plane thermopower and power fac- 
tor are presented at a electron/hole charge carrier con- 
centration of N = 3 x 10 19 cm~ 3 . The p-type ther- 
mopower shows only moderate dependencies on the SL 
period at all temperatures with the anisotropy slightly 
favouring the cross-plane part S ± . Under electron dop- 
ing the dependence of the thermopower on the SL pe- 
riod is more pronounced, which can to some amount 
be assigned to the quantum well effects which occur in 
the conduction band. At higher temperatures the ther- 
mopowers anisotropy is distinct larger. The latter was 
shown before for the bulk materials 18 . Due to the fact 
of the thermopower, as well as the electrical conductiv- 
ity being clearly smaller under electron doping than hole 
doping (cf. Figure 3) we find the largest values of PF 
for the p-type SLs. Here the largest was found to be 
PF y = 42 /iW/cmK 2 at 500 K for (Bi 2 Te 3 ) a; /(Sb2Te3)i_ :E 
SL with x = ye. For bulk Bi2Tc3 and Sb2Te3 we state 
maximum PF = 44/iW/cmK 2 and 49/iW/cmK 2 at 
400 K and 550 K, respectively. Due to previously dis- 
cussed conductivity anisotropy a u/cr±, the cross-plane PF 
is strongly suppressed and the maxima arc shifted to 
higher temperatures. 



C. Towards figure of merit 

With the electronic transport properties discussed in 
the previous sections, we are now going to focus on the 
electronic and lattice part contribution to the thermal 
conductivity n e i + K p h to give some estimations on the 
FOM. As has been stated before, the main benefit from 
a superlattice structure for the FOM is expected from 
a reduction of the cross-plane thermal conductivity at 
retained electronic transport properties. Today, the re- 
duction of the cross-plane lattice thermal conductivity in 
thermoelectric supcrlatticcs has been widely and success- 
fully proven 72-75 . 

In the past thermal conductivity reduction in crys- 
talline or polycrystalline bulk thermoelectric materials 
was traditionally achieved by alloying. However, one 
reaches the so-called alloy limit of thermal conductivity, 
which has been difficult to surpass by nanostructuring 2 . 

Nevertheless, for Bi 2 Te3/Sb2Te3 SLs cross-plane lat- 
tice thermal conductivites of k p h — 0.22 W/mK were re- 
ported for certain SL periods, which is a factor of two 
below the alloy limit 8 . It is obvious, that at thermoelec- 
tric relevant charge carrier concentrations and tempera- 
ture ranges, the electronic contribution n e i can be in the 
same order of magnitude. 

Therefore Figure 7 shows the room temperature dop- 
ing dependent electronic part of the thermal conductiv- 
ity, in the in-plane (thick dashed lines, right scale) and 
cross-plane direction (thin dashed lines, right scale), for 
bulk Bi 2 Tc 3 , to give insight in the principle dependen- 
cies. Furthermore, the Lorcnz function defined via Eqs. 3 
and 4 as Lj. ii = Keix,\\ ■ ( a ±.\\ ' T) -1 is shown for the 
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FIG. 6. (color online) Temperature dependence of thermo- 
electric transport properties for (Bi2Te3);r/(Sb2Te3)i_a; su- 
perlattices. Shown are (a) in-plane thermopower S» and 
(b) cross-plane thermopower S ± , as well as the correspond- 
ing power factors (c) PF|| in the in-plane and (d) PF^ in 
the cross-plane direction. The charge carrier concentration is 
fixed to N — 3 x 10 19 cm -3 and different line types correspond 
to different superlattice periods. Blue lines refer to electron 
doping, while red lines refer to hole doping. 



in-plane (thick solid line, left scale and color code) and 
cross-plane part (thin solid line, left scale), respectively. 
As can be seen, K e i minimizes for energies near the band 
edges. Here, at N w 3 x 10 18 cm -3 , the thermopower S 
maximizes at appropriate values for the electrical conduc- 
tivity a, hence the second term in Eq. 4 increases leading 
to small values for « e ;. At small intrinsic charge car- 
rier concentrations, the chemical potential shifts into the 
gap and the total thermopower is strongly reduced due 
to bipolar diffusion. This leads to an enhanced contri- 
bution to the electrical thermal conductivity at intrinsic 
charge carrier concentrations and is known as the bipolar 
thermodiffusion effect 76-78 . At charge carrier concentra- 
tions of N — 3 x 10 19 cm -3 we find K e i\\ to be about 
0.6 — 0.8 W/mK for n/p-type bulk Bi 2 Te 3 in very good 
agreement with experimental (cf. green, open circles in 
Fig. 7) and theoretical results 17,23,79 . The cross-plane 
component of n e i is substantially smaller, especially for 
n-typc conduction, reflecting here the electrical conduc- 
tivity anisotropy discussed earlier. The bipolar thcrmod- 
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FIG. 7. (color online) Lorenz function L (solid lines, ref. to 
left scale) and electronic contribution K e i to the total thermal 
conductivity (dashed lines, ref. to the right scale) in depen- 
dence on position of the chemical potential /i for bulk Bi2Te3 
for the in-plane (thick lines) and cross-plane (thin lines) trans- 
port direction. The Lorenz function is related to the metallic 
limit L = 2.44 x 10 _s Wn/K 2 . Plotted on to the Lorenz 
function in the in-plane direction is a color code referring to 
the charge carrier concentration given by Eq. 5. The red cross 
emphasizes the change from n to p doping, respectively. The 
temperature was fixed to 300 K. Thin vertical dash-dotted 
lines emphasize the position of the chemical potential for a 
charge carrier concentration of N — 3 x 10 19 cm~ 3 under n 
and p doping (blue and red line, respectively). The CBM 
is located at 0.105 eV. Green open circles show experimental 
values from Ref. 76 for K e i ii for an n-type Bi2Te3 single crystal. 



illusion effect is furthermore responsible for the suppres- 
sion of the Lorenz function to values below the metallic 
limit L Q (L = 2.44 x 10" 8 Wfi/K 2 ) for values of the 
chemical potential near the band edges (cf. Figure 7 
solid lines, right scale). At optimal charge carrier con- 
centrations of N = 3 x 10 19 cm~ 3 L|| « 0.7Lo under hole 
doping (red dashed dotted lines) and Lii « O.8L0 under 
electron doping (blue dashed dotted lines) can be found. 
For the cross-plane Lorenz function L± w L is stated 
at the same amount of n/p-type doping. Reaching the 
intrinsic doping regime the Lorenz function reaches sub- 
stantially large values of Lu ~ 6.5Lo and L± ks 8Lq. Such 
a behaviour has been described in literature 17 ' 80 and can 
have consequences for the determination of the thermal 
conductivity. The Lorenz factor is generally used to sep- 
arate K e i and K p h- At thermoelectric advisable charge 
carrier concentrations applying the metallic value Lo to 
determine the lattice thermal conductivity could lead to 
an overestimation of the electronic thermal conductiv- 
ity, and consequently to underestimation of the lattice 
contribution. The Lorenz function of thermoelectric het- 
crostructurcs can show further anomalies which arc dis- 



cussed in detail in a forthcoming publication 81 . 

Experimental findings for the lattice part n p h of 
the thermal conductivity are added to the calculated 
electronic contribution n e i to present some estima- 
tions on the cross-plane FOM ZT±. In particu- 
lar n p h.j_ = 1.05 W/mK, K p h,j_ = 0.96 W/mK and 
Kph.± = 0.22 W/mK at room temperature were used 
for bulk Bi 2 Te 3 , Sb 2 Te 3 and the Bi 2 Te3/Sb 2 Te3 SLs 8 , 
respectively 82 . 

Recently Winkler et al. 26 measured for a p-type 
(Bi . 2 Sbo.8)2Te 3 /Sb 2 Te3 SL the total cross-plane thermal 
conductivity n± to be about 0.45 — 0.65 W/mK for dif- 
ferent annealing temperatures. The values were obtained 
within a time-domain thermal reflectance (TDTR) mea- 
surement and are in very good agreement to our calcu- 
lations, which are displayed in Figure 8(c). Compared 
to the original experiments by Vcnkatasubramanian et 
al. 12 , the values of the total cross-plane thermal conduc- 
tivity are smaller. This stems to a large extent from 
the fact that a strong electrical conductivity anisotropy 
a \\/ax is apparent and not vanishing to cr u/cr_ L ~ 1 as 
proposed in Ref. 12. Hence not only a ± but also K e i ± 
is noticeably suppressed. Furthermore in Refs. 8 and 
12 it was suggested, that mirror- like SL interfaces lead 
to potential reflection effects and reduce n p h.± very ef- 
ficiently. Beside that, Touzelbaev et al. 24 showed that 
an existing interface roughness will additionally decrease 
n p h,±- Beside the nano-crystallinity of the samples, such 
an additional interface roughness is most likely provided 
by intcrdiffusion effects at the interfaces introduced by 
the growth of Bi 2 Te 3 /Sb 2 Tc3 SLs within the concept of 
"nano- alloying" 2 ^ 26 . 

In Figure 8 the room-temperature cross-plane prop- 
erties of the total thermal conductivity and the related 
cross-plane part of the FOM are shown, for electron 
(a),(b) and hole doping (c),(d), respectively. As could 
be expected the strong quantum well effects in the con- 
duction band lead to quite small values for ZTj_ < 0.2 
under electron doping. Here a reduction of the SL total 
thermal conductivity of about a factor 5-6 compared to 
bulk is impeded by an accompanied reduction of PF^ by 
a factor of 10-20 (cf. Figure 4(c)). Thus no benefit for 
the FOM in the cross-plane direction can be revealed un- 
der electron doping. At hole doping the situation is more 
advantageous. Beside the total thermal conductivity (cf. 
Figure 8(c) and Figure 7) being somewhat larger than 
under electron doping, the electronic transport proper- 
ties remain bulk-like for all SL periods (cf. Figure 4(f)) 
and hence an enhancement for ZTj_ by a factor of 2-3 
can be achieved under hole doping of N = 3 x 10 19 cm" 3 . 
We state values ZT ± « 0.7 - 0.9 for different SL peri- 
ods, while an higher amount of Sb 2 Tc3 in the SL leads to 
larger values for ZT, despite no distinct influence of the 
SL period on the FOM could be found. As a supplement 
experimental data 64 for k± and ZTj_ of the mixed single 
crystal series (Bi x Sbi__a;) 2 Te3 is shown as green down- 
ward triangles in Figure 8(c) and (d). The hole doping 
varied steadily between N — 1 x 10 19 cm -3 for Bi 2 Te3 to 
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FIG. 8. (color online) Absolute values of cross-plane thermoelectric transport properties for (Bi2Te3) a: /(Sb2Te3)i_ a ; superlattices 
in dependence on the superlattice period. Shown are (a),(c) the total thermal conductivity n e i + K p h and (b),(d) the cross- 
plane figure of merit ZTj_. The temperature is fixed to 300 K and results for three different charge carrier concentrations are 
compared, (a), (b) refer to electron doping, while (c), (d) refer to hole doping. The electronic part n e i was calculated, while 
the lattice part K p h was taken from literature 8,12 . As a reference point experimental results for the mixed single crystal series 
(Bi^Sbi-z^Tes at comparable material composition are shown as green downward triangles 64 . 
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FIG. 9. (color online) Temperature dependence of 
the thermal conductivity and figure of merit for the 
(Bi2Te3) a; /(Sb2Te3)i_ a ; SLs. Shown are (a) total cross-plane 
thermal conductivity n± = « e i,± + n- P h,± and (b) cross-plane 
figure of merit. The charge carrier concentration is fixed 
to N = 3 x 10 19 cm -3 and different line types correspond 
to different SL periods. Blue lines refer to electron doping, 
while red lines refer to hole doping. The electronic part K e i 
was calculated, while the lattice part k p k was taken from 
literature 8,12 . 



N = 9 x 10 19 cm -3 for Sb2Te 3 . One easily abstracts that 
ZT± values of the optimal doped SLs and the mixed sin- 
gle crystal series show clear similarities considering the 
dependence on the materials composition, as well as the 
the absolute values of ZTj_. 

To extend our findings at room-temperature, in Figure 
9(a) and (b) temperature dependent results for the total 



cross-plane thermal conductivity and cross-plane FOM 
are shown for the (Bi 2 Te3) 2: /(Sb2Te3)i_ a ; SLs at elec- 
tron/hole concentration of N = 3 x 10 19 cm -3 . Within, a 
conventional 1 /t dependence for the lattice thermal con- 
ductivity was assumed 83 , while the calculated electronic 
part K e i is temperature and doping dependent, per se. 
We note that for thermoelectric SLs no clear tendency 
on the temperature dependence of n p h can be revealed. 
However, conventional 1 /t dependence, as well as temper- 
ature independent n p h were found experimentally 72 ' 73 ' 84 . 
Models show, that K p h should diminish at low periods 85 , 
while experiments reveal a saturation towards the alloy 
limit for SL periods below 50 A 8 . 

As can be seen from Figure 9(a) k± takes a minimum 
at about 300 - 400 K for the Bi 2 Te3/Sb 2 Te3 SLs. This 
behaviour is dictated by the electronic contribution to 
kj_ and supported by the 1 /t dependence of the lattice 
part of k. At low temperature the chemical potential 
is located in the bands an thus a moderate contribu- 
tion to n e i,± is obtained. With increasing temperature 
K e i,± slightly decreases as the chemical potential shifts 
towards the band edges, then reaching minimal K e /,j_ for 
chemical potential positions at the band edges (cf. Fig- 
ure 7). At elevated temperatures the bipolar contribu- 
tion leads to an enhanced electronic contribution to 
which then clearly dominates the 1 /t dependence of the 
lattice part of n leading to large values of the total ther- 
mal conductivity. The influence of the electronic con- 
tribution is more pronounced in the SLs compared to 
bulk, as here K. e i,±. > k p h,±- Combining these results 
with the temperature dependent power factor PF^ dis- 
cussed in Figure 5(d) we find the temperature depen- 
dence on the cross-plane FOM as presented in Figure 
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9(b). Concentrating on the more promising p-typc SLs 
we state maximized values for ZT_l clearly above unity 
for temperatures of 400 — 500 K. The largest cross-plane 
FOM is found to be ZT ± = 1.27 at about 470 K for a 
(Bi 2 Te3) 2; /(Sb2Te3)i_ a; SL at x = 1 /e and a hole concen- 
tration of N = 3 x 10 19 cm~ 3 . We want to mention that 
in the experiments of Ref. 12 the maximal ZT± under 
hole doping was stated at a SL period of x = l /6, too. 
The best value for an n-type SL is ZTj_ = 0.25 at about 
360 K at a SL period of x = 5 /6. 

V. CONCLUSION 

The anisotropic thermoelectric transport properties of 
Bi2Tc3/Sb2Tc3 supcrlatticcs at different superlattice pe- 
riods is presented to get insight into the physical mech- 
anisms which are responsible for the path-breaking ex- 
perimental results with ZTj_ = 2.4 at room temperature 
obtained by Ref. 12. Several aspects added up to obtain 
those very high ZT values in experiment, (i) In-plane val- 
ues of the electrical conductivity and thermopower were 
found to be larger than in the bulk systems at compara- 
ble charge carrier concentrations, (ii) An elimination of 
the electrical conductivity anisotropy < r \\/<r±, which is ap- 
parent in both bulk systems, was found in supcrlatticcs 
at certain periods. (Hi) The lattice part of the thermal 
conductivity was reduced below the alloy limit due to 
phonon-blocking at the superlattice interfaces 8,24 . 

Even though taking into account the most optimistic 
value K p h.± = 0.22 W/mK for the lattice part of the ther- 
mal conductivity, we found the cross-plane figure of merit 
for the best p-type superlattice to be ZTj_ = 0.9 at room 
temperature and slightly enhanced to ZT^ = 1.27 at ele- 
vated temperature. However, this is more than a factor of 



two worse than experimentally revealed and is caused by 
the fact, that within the presented ab initio calculations 
the findings (i) and (ii) could not be confirmed. We want 
to add, that in first consistent experiments by Venkata- 
subramanian et al. 62 a room-temperature ZT^ = 1.2 
was proposed for a non-symmetrical Bi2Tc3/Sb2Te3 su- 
perlattice with a period of 30A. 

For the in-plane transport properties of S, a and PF 
we can state values comparable to bulk for the p-type 
and n-type superlattices, which is in agreement to recent 
experiments 21,26 . Furthermore for the p-type superlat- 
tices a conservation of the bulk transport anisotropics 
is found, but in no case a reduction, while under elec- 
tron doping strong quantum well effects due to con- 
duction band offsets lead to large transport anisotropics 
a \\/<r± > 10 and suppress the cross-plane thermoelectric 
transport notably. 

Concluding, the experimentally found remarkable ther- 
moelectric transport properties in Bi2Te3/Sb2Te3 super- 
lattices could not be revealed by detailed band structure 
effects. An ongoing issue will be to clarify, whether scat- 
tering effects caused by defects and lattice imperfections 
could give a considerable leap forward to understand the 
enhanced thermoelectric efficiency in Bi2Tc3/Sb2Te3 su- 
perlattices. 
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